Matter-field theory of the Casimir force 
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A matter-field theory of the Casimir force is formulated in which the electromagnetic field and 
collective modes of dielectric media are treated on an equal footing. In our theory, the Casimir 
force is attributed to zero-point energies of the combined matter-field modes. We analyze why some 
of the existing theories favor the interpretation of the Casimir force as originating from zero-point 
energies of the electromagnetic field and others from those of the matter. 

PACS numbers: 12.20.Ds, 42.50.Ct, 42.50.Lc 



I. INTRODUCTION 



P-f 

, It is well known that two conducting neutral plates placed at a small distance attract each other. This phenomenon 
' is often called the Casimir effect because of the celebrated formula he derived for the force between two perfectly 
conducting plates 0] . Casimir calculated the sum of the quantum-mechanical zero-point energies of the normal modes 
of the electromagnetic (EM) field between the two metal plates, and showed that the total energy depends on the 
distance between the plates. The spatial derivative of the sum gives what we now call the Casimir force. The 
generalization of the Casimir effect to the case of two dielectric plates was made by Lifshitz ^ , and his formula was 
rederived by van Kampen, Nijboer, and Schram j^^] using a technique attributed to Casimir. The Casimir force is 
essentially a long-range van der Waals force, where one cannot ignore the delay caused by the finiteness of the velocity 
of light. This retardation effect is not always a small correction and may alter the very nature of the force. For 
example, the Casimir force between two perfectly conducting hemispherical shells has been shown to be repulsive [0 . 

Experimental studies of the Casimir force started long ago and have culminated in a recent high-resolution 

measurement by Lamoreaux B. A closely related phenomenon, the retarded attraction force between an atom 



and a metal plate (the Casimir-Polder force |10| ) , was also measured successfully by Sukenik et al. 1 1 1 . With further 
improvements, it will hopefully be possible to confirm experimentally the theoretical predictions of several corrections, 
such as finite-temperature effects |l^,|ll||l4| and radiative corrections ]l5| . 

The scope of the research of the Casimir effect covers many areas of physics and other fields, ranging from biol- 
ogy to cosmology and elementary particle physics. Reviews and extensive references are available, for example, in 
Refs. 



The Casimir effect has usually been attributed to zero-point fluctuations of the EM field. However, Schwinger and 



his collaborators 14 2C|] showed that the Casimir effect can be derived in terms of Schwinger's source theory, which has 
no explicit reference to vacuum-field fluctuations of the EM field. Milonni and Shih have recently developed a source 
theory of the Casimir effect within the framework of conventional quantum electrodynamics pi| . In this theory, the 
Casimir force originates from quantum fluctuations of atomic dipoles in the dielectric, and the EM field plays only 
?3 i the passive role of mediating interactions between those dipoles. 

It seems like it is only a matter of taste whether we attribute the Casimir effect to the quantum nature of the 
EM field or to that of the matter. Accordingly, the following questions naturally arise: Do the field and the matter 
really play (or to what extent do they play) symmetrical roles in this problem? Why do some approaches emphasize 
the quantum nature of the EM field and others stress that of the matter? To answer these questions, it would be of 
interest to discuss this problem from a standpoint that has no preference for the field or for the matter. Unfortunately, 
existing theories do not suit this purpose because they all invoke in one way or another the Maxwell equations in 
which the degrees of freedom of the matter are in advance embedded in the frequency dependence of the dielectric 
response function. 

The primary purpose of this paper is to propose a matter-field theory of the Casimir force in which the matter and 
the field are treated on an equal footing. Our strategy is to explicitly diagonalize the matter-field Hamiltonian which 
is quadratic in its dynamical variables. Because relevant physical quantities are then expressed in terms of eigenvalues 
and eigenvectors of the full Hamiltonian, all physical effects allow unambiguous interpretation. 



*To be published in Phys. Rev. A 
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This paper is organized as follows. Section || shows how the Casimir force is derived, starting from the Lagrangian 
describing the interaction between the EM field and collective modes of the matter. Section [II proves that this 
derivation leads to the same results as existing theories [[l],p3JI|,|l4|,H(ni . It will be shown that in our treatment both 



field and matter contributes zero-point energies. Section pV|analyzes the previous derivations of the Casimir force 
and points out that an inherent asymmetry indeed exists in the problem between the field and the matter, which 
favors the interpretation of the Casimir force as originating from zero-point fluctuations of the EM field in some 
formulations ]3|,|| , and from those of the matter in others pi| . Section ^ summarizes the main conclusions of this 
paper and describes our answers to the above questions. 



II. MATTER-FIELD THEORY 

A. Formulation of the problem 

Consider a system in which two dielectric slabs are separated by the vacuum of the EM field as shown in Fig. [I]. We 
assume that the entire system V is enclosed by a perfectly conducting metal, and denote the boundary of the system 
on the metal as S. The region V m occupied by the dielectric is surrounded by S and by the interface S m between 
the dielectric and the vacuum. The following discussion does not depend on concrete shapes of V and V m , but we 
require that any closed loop or any closed surface in V can continuously shrink to a point without going outside of V. 
The condition concerning closed loops ensures that any irrotational vector function Q(r) [i.e., V x Q(r) = 0] has a 
scalar potential (j>(r) with — V0(r) = Q(r). The condition concerning closed surfaces ensures that if Q(r) has a scalar 
potential and Q _L S, Q(r) has a scalar potential <p that vanishes everywhere on S. We impose the same requirements 
on the topology of V m . 
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FIG. 1. Two dielectric slabs (region V m ) are separated by the vacuum. The entire system V is enclosed by a perfectly 
conducting metal. The boundary of the system on the metal is denoted by S, and the boundary between the dielectric slabs 
and the vacuum of the EM held is denoted by S m - 

In a field theory of the Casimir force Q|, normal- mode frequencies of the "EM field" are determined from the 
Maxwell equations with proper boundary conditions, and the degrees of freedom of the dielectric slabs do not appear 
explicitly because the effect of the presence of the dielectric is incorporated into the theory through the dielectric 
response function e(u>). In a matter theory of the Casimir force pl| , the Green function of the EM field, which 
mediates the interaction between atomic dipoles, is determined from the same dielectric response function. Because 
our goal is to treat the EM field and the matter on an equal footing, we first consider the EM field and the dielectric 
media separately, and then introduce the interaction between them. 

The Lagrangian Lgeid of the free EM field is written as 

L 'W/"(?' E l 2 -5> 2 )' W 

where we assume that the EM field is only subject to a perfectly conducting boundary S. 

We take the following simple model for the dielectric: It consists of various kinds of particles whose species are 
labeled by j with number density rij. The jth species has mass rrij and charge ej, and particles belonging to this 
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species are bounded around their equilibrium positions by a restoring force characterized by frequency u)j. In order 
to describe their collective motions, we use a complete orthonormal set {/3j(r)} of mode functions in the region V m 
that is occupied by the dielectric. They satisfy 



(2) 



where is the Kronecker's delta. We take the boundary condition on S as Pi || S and that on S m as /3, || S m . 
Note that imposing a boundary condition does not alter the physics in our model because we assume no interaction 
(except that mediated by the EM field) between the particles. Under this boundary condition, the set {Pi] is uniquely 
classified into the set of transverse vector fields {p\ '} satisfying V • p\ =0, and the set of longitudinal ones {P\ } 
satisfying V x Pi' = 0. 



The polarization density Pj (r) of the j th species can be expanded in terms of these bases as 



P i( r) = ^£/3P(r)xf >(r)X<f>, 



(3) 



(T) (L) 

where XL and are the generalized coordinates that describe the collective displacements of the jth species. 

The corresponding current density jj (r) and the charge density pj (r) are given by 

j,(r) = P,(r) = ^T J2PP(r)Xp + P[ L \r)X 

i i 

and 

Pj (r) = V ■ P,(r) = -^T J2 V • P\ L) (r)X^ 



(L) 
ij 



(4) 



(5) 



In Eq. (||), we have neglected nonlinear terms in X^J' L \ which give rise to the term v x B of the Lorentz force. They 
may provide a relativistic correction to the Casimir force, which, however, does not arise in the lowest nontrivial order 
we are concerned with. The total kinetic energy of the jth particle is written as follows: 



T; 



dv—n-im-i 
2 3 3 



■'3"'] 



ij I 



(6) 



where 



j i>j 



= J. 



e j n 3 / e o m j is the plasma frequency. The Lagrangian for the dielectric L D 



is thus given by 



<T)2 

ij ^PJ 



,, 2 Y^ 2 



,(L)2 

2uj 2 - y '^ 3 

ij PJ 



,,2 v(L)2} 



The total Lagrangian L for the whole system is the sum of Lfi e id, -^matter, and the interaction part: 

L = Lfloid + ^matter - V" / dv [pj(r)cp(r) - ]j(v) ■ A(r)] , 

Jv m 



(7) 



(8) 



where ip(r) and A(r) are respectively the scalar and the vector potential of the EM field. The corresponding Hamil- 
tonian in the Coulomb gauge is derived as follows: 



ij 



dL 



dX. 
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where E^ T ^ and E' l ) are respectively the transverse and the longitudinal parts of the electric field; -Py = dL/dX^ 

and P}^ = dL/dX^ arc the generalized momenta conjugate to X^p and X^\ respectively. 

The transverse part of the EM field can be decomposed into the normal modes in the usual manner. Consider a 
complete orthonormal set of transverse mode functions {ol\ (r)} in V that satisfy 

A 2 af)(r) =-%P(rl (10) 



J a f\r)- a ^(r)dv = S tJ . (11) 



(T) (T) 

Because S is the boundary on the perfect conductor, al should satisfy the boundary condition a\ _L S. The 
transverse fields A(r) and E^ T ^(r) can be expanded with this base as 

A(r) = JL^a( T )(r) gi (12) 



and 



E( T )(r) = 4E« (T) WP.. (13) 



'A 

where qi and pi are canonically conjugate pairs. The transverse EM field part of the Hamiltonian ^ then reduces to 
a collection of harmonic oscillators: 



£*(f | " m|>+ 5> , )-sG , * + 5 n 



hi ] ■ I") 



In the Coulomb gauge, the longitudinal part E^ L ^(r) of the electric field satisfies the Poisson equation eoV-E^ = p, 
where the total charge p is given from Eq. (||) as 

p = E Pi = -v^E v • $t\*) x f ( 15 ) 

i y 

in V m and zero elsewhere. Let us define the longitudinal mode function a\ (r) in V that satisfies the boundary 
condition a ^ _L S and 

V-^ ) (r) = (7-^ ) W in >. (16) 
[0 otherwise. 

The solution to Eq. (|l^) is given as 

f3f\r) -E4 T) W / aL T) (r')-#V)*/ in 7 m , 



"f } (r) = < 



(17) 

-Vaf } (r) / a { p{r')-(3[ L \r')dv' otherwise. 



Note that oq (r) is neither normalized nor orthogonal. With this base and x\^\ E^ L )(r) can now be expanded as 
follows: 

EW(r) = -^5>f)(r)Xf\ (18) 

v u y 

The Coulomb-energy part of the Hamiltonian (0) is 



cx^(v)- a [ L \v)dv. (19) 



4 



Using Eq. (|T^) and the orthogonality between ol\ (r) and or k ' (r) , the integral on the right-hand side of Eq. ([H 
can be rewritten as follows: 



Combining Eqs. 



/ af ) (r).4 i) (r)dw= / ^(r) ■ a [ L) (r)dv. 



(20) 



fll2| ), (|14|), (|19|), and (|20|), we obtain the following expression for the total Hamiltonian: 



Oj Y (T)2 



2a; 2 y 



E 



0j _ X (L}2 



2lo 2 ■ 13 
pj 



/3^(r).«^(r)d W . 



(21) 



B. Plasmons 

In this section, we discuss the physical implications of the longitudinal displacement represented by {p[ L \r)}. Let 
us consider another set of complete orthonormal modes {/3j(r)} in V m that satisfy different boundary conditions on the 
boundaries, namely, f3 i _L 5 and f3 t _L S rn . Under these new boundary conditions, the set is uniquely separated 

into the set of transverse vector fields {f3 i } satisfying V • f3 i = 0, and the set of longitudinal ones {(3 i }, which 
have scalar potentials </>i(r) (p i = — V0j) that vanish on S and S m . 
The orthogonality between /3 fc and /3^ T ' ) is shown as follows: 

^ /3 4 (T) (r) • (r)dv = - j n • /3f ' (r)& (r)dS = 0, (22) 

where the surface integral runs over the surface of V m , and n is the unit vector normal to that surface. In contrast, 
@\ "* and (3^ are not necessarily orthogonal. This implies that a vector longitudinal under the boundary conditions 
(3i || S, S m is not necessarily longitudinal under the boundary conditions /3, _L S, S m . Let us choose the base 
as a union of two such orthogonal sets, {(3^} and {/3^}, that f3^ can be expanded with {f}\ }, and with 
{/3^ }. Associated with (3^ and /3^ s \ we may introduce ct^ and ot[ s ^ through the relation (|l7|). 

Since /3j; has a scalar potential <pf' ((3^ — — V0^) that vanishes on the boundaries 5 and £? m , part of the 
integral ( |20| ) that includes the mode can be further simplified to be 

/ f3 { l\v)-a.[ L \v)dv= [ ^f ) (r)V-4 i) (r)d U = / ^ 6) (r)V • (3 { k L) (r)dv 
JV m JV„ Jv m 

(3f\v)-(3[ L \v)dv, (23) 

where we used Eq. (|l|) in deriving the second equality. Because of the orthogonality condition between pf 1 and 
we obtain 

/ (3?\v)-ct ( k s) (r)dv = (24) 
Jv m 

and because of the orthonormality of {/3^ }, we obtain 
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P\ h \r).cx^(r)dv = 6 ik . (25) 



From these, the static Coulomb-interaction part of the Hamiltonian can be decomposed into two parts: 
^(e^) (e*#) X m ^(r) a«(r) dv + ii: ^) 



(26) 



This result shows that the Coulomb interaction between the modes xjj does not depend on the boundaries (bulk 
plasmons). Further, these modes do not interact with the transverse EM field since 

J af ) (r) • /3< b) (r)dv = - J n ■ aj T) (r)^ 6) (r)dS = 0, (27) 

where the surface integral runs over the surface of V m . The bulk plasmons thus do not contribute to the Casimir 
force. 

(s) 

The modes X>, describe the surface plasmons whose resonance frequencies do depend on the parameters of the 



surface of V m and its outside, through the integral in Eq. (E6J). The surface plasmons are also affected by the 
surroundings through the interaction with the transverse EM field, for their scalar potentials by definition should not 
vanish everywhere on the surface of V m and the integral J v c*^ (3% dv does not vanish unlike Eq. (p7|). 



C. Combined oscillations 



After separating the longitudinal modes into the two types of plasmons, the total Hamiltonian is now written as 
follows: 



where 



dv 



^Oj y(T)2 

2lu 2 - ij 

P3 



fOj y(s)2 

2uj 2 ■ lj 
pj 



EG 

i x 

E A 



2 11 



pipiw 

2 « 



"0j y(b)2 



2uL ij 



K, 



k = I ^(r) ■ a[ s '(r)dv 



(3\ s) (r)-a [ k s) (r)dv. 



(28) 



(29) 



For simplicity, we introduce new conjugate variables qi and pi by canonical transformation pi = —qi and qi = pi 



If we define column vectors x =' {{X\* J }, {X}*'}, {&}, {X% } }) and p =' ({P-- '}, {P-j }, {pi}, {P}j }), the total 
Hamiltonian can be written in a compact form as follows: 



H = i*pTp + -*xFx, 



(30) 



G 



where T and V are symmetric matrices. Since T is positive definite, it has a decomposition T = A 2 by a symmetric 
regular matrix A. The positive semidefinite Hermite matrix AVA can be diagonalized by an orthogonal matrix 
U as t UAVAU = D, where D is a diagonal matrix with no negative elements. A transformation to new bases 
X = t (Xi,X2, ■ ■ .) and P = t (Pi,P2, ■ ■ ■) defined as X = t UA~ 1 x and p = A -1 UP is a canonical transformation 
since the commutation relations are preserved: 

[XhPj] = [Y^CUA^Uxk^UA)^] = ^UA-^ikCUA^ix^pt] 

k I kl 

= ihJ^CUA-^kiAU)^ = iMij. (31) 

k 

This transformation diagonalizes the Hamiltonian (|3^) into a sum of independent harmonic oscillators: 

H = i*PP + i*XDX = £ \{P 2 + utXf) + 2 \{Pi b)2 + (32) 

i k 

where we have separated the contribution of bulk plasmons and suffixed their variables with superscript (b) since they 
are decoupled from the others. 

Now we second quantize Eq. (^) to obtain 

H = £ huiMU + ~) + E ^i 6) (4 b)t 4 6) + 5), (33) 

where we introduced creation and annihilation operators in the usual way. The energy of the ground state (H)q is 

(H) = E Q + E^ ^X>« + ^E ^ ( 34 ) 

i k 

The Casimir force is derived from the change in the ground-state energy with respect to an infinitesimal displacement 
of the dielectric slabs. As we have seen, the energy contributed by the bulk plasmons is independent of this 
displacement. Thus the Casimir force originates from Eq — the sum of zero-point energies of the harmonic-oscillator 
modes X;. The relation X = t UA^ 1 x. implies that Xi is a linear combination of g,, xj-j, and X^ . This means that 
Xi represents a combined mode of the EM field and the collective modes of charges in the dielectric. Therefore, in our 
approach, the Casimir force is attributable neither to the change in zero-point energies of the genuine EM field nor 
to that in zero-point energies of the genuine matter, but to that in zero-point energies of the combined matter-field 
modes. 



III. EQUIVALENCE TO MAXWELL EQUATIONS 



In this section, we show that the diagonalization of the Hamiltonian (g8|) amounts to solving the Maxwell equations 
associated with a proper dielectric response function e(u>) under appropriate boundary conditions. 

Consider the equations of motion derived from the Hamiltonian (p8|). After eliminating % P^J\ p}j \ an d Pij\ 
these equations read 

Pi 



^+E</ ^-/C^ + E^?/ ^-^dv, (35) 

n? = -<4? - Eft / ^ • «r ^ (36) 

x S=-^ x S-4jj:4% (38) 
3' 

The diagonalization of the Hamiltonian is equivalent to finding the solution to the above set of equations that oscillates 

at frequency u. For such a solution, the second-order time derivative may be replaced by — uj 2 . Equation (|3^) then 
reduces to 
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Pi ~ 2 02 



where we used the notation 



4 M) -E4f ] (M=T,s,b). 



From the remaining equations (|36|)-(|38|), we obtain after summation over the suffix j 



XW = (1- C ( W )) 



= (l-e(u;))xf, 

where we have introduced the dielectric response function e(u;) defined as 

9 



- Ufij ' 



Let us introduce the electric field E(r) defined as 

E(r) = fcftoTW + E*«a«(r) + E^^W) 

v e ° V i k I J 

Then it is easy to show that Eqs. j4]])— (|43|) have the same form as follows: 

4 M) = >/^(c(w) - 1) / P[ M) • Vdv (M = T, s, 6). 



Substituting this into Eq. (|39|) and noting the completeness relation 

i ft / 

where 1 denotes the identity matrix, we obtain 

Substituting Eqs. (^) and (^|) into Eq. ( [45| ) yields the self-consistent condition for the electric field E(r): 
E(r) = [1 - e(w)] (£ _J^ a < T >( r ) ^ af V) • E(r')<fo' 

+ E«i S) (r)/ y /3i s) (r')-E(r')^' + E«| b) w/ y $ V) • E^dt^ . 

This equation can be simplified if we convert it into a pair of equations as follows. First, by taking the diver 
both sides we obtain 



V-E(r) = [l-e(o;)]f£V-a«(r) J v /3iV)-E(r')^' 



Substituting Eq. (|16|) into Eq. (|50|), we obtain 



V-E(r) = [l-e(«)]V. ( £ (r) jf ^ • E(r')^' 
+ £/3f ) (r) jf /3[ b) (r')-E(r')^9(r), 



(51) 



where O(r) is the unit step function defined by 



®( r ) | otherwise. 

Using the completeness relation (^) and noting that /3- T ^(r)9(r) is a transverse function in V, Eq. fl5l| ) is simplified, 
giving 

V.E(r) = [l-e(w)]V-[E(r)e(r)]. (53) 
On the other hand, by operating (uj 2 — c 2 rotrot) on both sides of Eq. (ffi^ ) and using Eq. (px|) , we have 

(w 2 - c 2 rotrot)E(r) = [1 - e(w)]w 2 ^ a| T) (r) j af ) (r') ■ E(r')^' 

+ J24 S \r) [ /3i s) (r')-E(r')^' + E«| b) W / /fV) ' E(r')d«'V (54) 
Here, using Jl7j ) and ([l7|), we obtain 

£ «f > (r)«f } (r') + E «i S) (O + E «, W (*)ti b) ^) 

i k I 

= E e(r)/3«(r)/3<V) + £ 6(r)/3 ; (b) (r)/3 ; (b) (r') 

fe / 

+ E«f ) Waf ) (r')-E a r ) W / afVT^VWiV) 

= E 0(r)/3( s )(r)/3^(r') + £ 6(r)/3[ b) (r)/3 ; (b) (O 
fe / 

+ 5>f>(r) / cc?\r").(3£\r")dv"p£\r>) 

JVm. 



= 9(r)l<S 3 (r-r'), (55) 
where we have used the fact that Q(r")f3?"* (r") can be expanded with a^(r"). We thus obtain 

{uj 2 - c 2 rotrot)E(r) = [1 - e(tj)]u 2 E(r)e(r). (56) 
If we define the electric displacement D(r) as D(r) = eo{l + [c(ui) — l]0(r)}E(r), we find that Eqs. (53) and (p36|) 



are identical to the Maxwell equations, namely, 

V • D(r) = (57) 

and 

V x (V x E(r)) = — -D(r). (58) 

Our theory developed in Sec. |lj therefore yields the same sequence of eigenvalues u>i as Schram's method Q and 
eigenvalues of the bulk plasmon modes that are not relevant to the Casimir force. Our theory thus gives the identical 



result for the Casimir force: for example, when the interface S m consists of a pair of parallel planes with a distance 
d, the force per area is calculated to be 

f- poo />oo 



where s = \J j> 2 + ei£ — 1. 



' s +pe(i£) 



- 1 




a 2pt;d/c 



- 1 



(59) 



IV. FIELD THEORIES AND MATTER THEORIES OF THE CASIMIR FORCE 



In the preceding sections, we have developed a theory of the Casimir force by starting from the Lagrangian that 
describes both the EM field and the collective modes of charges in the matter. In this theory, the Casimir effect 
is attributed to a change in zero-point energies of the combined modes of the EM field and the matter. Thus, in 
our theory, both field and matter contribute to the Casimir force. The system considered here, however, has one 
asymmetry between the field and matter, which is obvious from Fig. [l]. The asymmetry originates from the property 
of the interface S m in Fig. [I]. It confines the matter in one side, but imposes no boundary condition on the EM field. 
It is difficult to imagine an opposite situation. That is, we can exclude the matter from some region, but we cannot 
clear away the EM vacuum. As will be discussed below, it is this asymmetry that makes some of the existing theories 
favor the interpretation of the Casimir force as arising from zero-point EM energies (field theories), and others favor 
the interpretation as due to zero-point fluctuations of the matter (matter theories). 

First, we consider a field theory of the Casimir force, by which we refer to the scheme of solving the Maxwell 
equations under appropriate boundary conditions. The relation between our theory and the field theory has already 



been discussed in Sec. [II. There we have shown that both theories actually calculate the zero-point energies of 
the same set of combined matter-field modes. Nevertheless these energies are sometimes incorrectly identified as 
zero-point energies of the genuine EM field, presumably because the normal modes can be specified by only invoking 
the Maxwell equations without any reference to the state of the matter. The degrees of freedom of the matter arc 
embedded in the dielectric response function from the very beginning in the field theory of the Casimir force and are 
therefore not manifest. In addition, in the simplest case of the two perfectly conducting plates, it would be easy to 
forget the matter because it occupies only the thin surface and the whole mode volume is filled with the pure EM 
field. As a matter of fact, however, within this thin surface region there exist rather complicated surface modes, as 
described in Sec. 0, which conspire to shield the EM field, producing the Casimir effect. In contrast, these normal 
modes can hardly be looked upon as genuine matter oscillations because it is difficult to construct a simple wave 
equation that contains only the matter variables when the system includes regions where no matter exists. 

A matter theory of the Casimir force, which is based on zero-point fluctuations of the matter, is formulated by 
Milonni and Shih |H] in terms of conventional quantum electrodynamics. Their scheme looks somewhat complicated, 
but the essential ingredient of the theory is the use of second-order perturbation theory in calculating the self-energy 
of linearly interacting harmonic oscillators. Here we will present the bare essentials of this theory by taking up the 
simplest example of a pair of interacting harmonic oscillators. 

The Hamiltonian for the system is written as 

H = H Q + V, (60) 

where 

Ho = Ku a (tfa + + fiuj b (b r b + , (61) 

V = v+atP +v-tfb + v* + ba + v*_tfa, (62) 

and the interaction V is assumed to be small. In this case, the first nonvanishing correction to the ground-state energy 
is of second order in V and is written as AE^ = (g^\V\g), where the state \g^) is the unperturbed ground state 
and \g) is the ground state, which is correct up to first order in V. Instead of using a time-independent perturbation 
theory, we can also obtain \g) by a time-dependent perturbation theory with the Hamiltonian Hq + Ve Jt , where the 
limit 7 —> +0 is implied. Hence 

|S) = {l- l -J° dt'e^'Ve^'e-^A \g^) 

= U\gW). (63) 
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Up to second order, we have 



AE® = (g i0) \VU\gM) = (g( 0) \U*V\gW) 

= \{g^\uWU\g^). (64) 

This corresponds to the interaction energy assumed in Ref. |2l]] . The following relations are useful for the following 
calculation: 

U^aU = a ^± — o f — — b, (65) 

0J a +U> b UJ a - LU b 

U r bU = b — — ^ — — a. (66) 

UJ a + LUb U b - UJ a 

Since the interaction V consists of the products of commuting operators, changing the ordering of operators in V 
should not affect the value of AE^ 1 . Let us take one particular ordering of 

V = v + cJP + v-dlb + v* + ba + vltfa, (67) 

and substitute Eqs. (|6^) and ( |66| ) directly into U^VU without changing the order of operators. In the resulting 
expression, the operators for mode a appear in normal order, and therefore do not contribute to the expectation value 
in the vacuum. Thus, AE^ in Eq. (J6J) is expressed only in terms of zero-point fluctuations of mode b as 

AE^ = ^+L( 5 (o)|SSt| 5 (o)). ( 68 ) 

UJ a + LU b 

If we take another choice of operator ordering for V, 

V = u+S t a + +v-a}b + v* + ab + v*_b^a, (69) 

the same procedure leads to 

A£(2) = ^+L/ ff (o)i aa tu(o)\. (70) 

LO a + UJ b 

In this case, the correction is attributable to quantum fluctuations of mode a. The second-order correction to the 
ground-state energy of a pair of linearly coupled harmonic oscillators can therefore be attributed solely to zero-point 
fluctuations of either one of the oscillators. 

Now, back to the original problem of the Casimir force. We notice [^lj that the force is derived from the change in 
the energy of the entire system when we add atoms in an infinitesimally thin layer next to the boundary S m |^| . The 
dominant contribution to this change arises from the interaction energy between the added atoms and the field modes, 
either of which may be described as a set of harmonic oscillators. Thus, as in the simple system discussed above, we 
can attribute the change in energy to the fluctuations of the added atoms alone, provided that an appropriate operator 
ordering is chosen (matter theory). We emphasize that the logical consistency of this interpretation hinges heavily 
on second-order perturbation theory as applied to a system of linearly coupled harmonic oscillators (linear-response 
theory). If effects of higher-order interactions are not negligible, it would be nontrivial to construct either genuine 
matter theory or field theory that is consistent with each other. 

If we choose another appropriate operator ordering, the Casimir force may be interpreted as arising from quantum 
fluctuations of the EM field modes, and we obtain a kind of field theory of the Casimir force. But the machinery is 
similar to the field theories described earlier, that is, the field modes are actually not the genuine EM field modes but 
the combined modes of the EM field and the matter. 

Second-order perturbation theory (or linear-response theory) thus allows the interpretation of the Casimir force as 
arising from quantum fluctuations of the genuine matter alone but not from those of the genuine EM field alone. The 
reason for this asymmetry lies again in the property of the boundary. That is, an infinitesimal displacement of the 
boundary only affects the mode functions of the genuine matter but has no direct effect on the mode functions of the 
genuine EM field. 
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V. SUMMARY 



The origin of the Casimir force has usually been attributed to zero-point fluctuations of the EM field, and sometimes 
to those of the matter. We have formulated the problem by starting from the Lagrangian that describes the EM field, 
collective modes of the matter, and their interaction. This approach makes it clear that the Casimir effect, in actual 
fact, arises from the change in zero-point energies of certain combined modes of the transverse EM field, the transverse 
motion of charges, and the surface plasmons. The zero-point energy of the entire system is therefore contributed by 
both the EM field and the matter. 

The Casimir force arises not from the zero-point energy per se, but from its change with respect to a virtual 
infinitesimal displacement of the plates. This change alters the mode functions of the matter alone, and docs not 
affect those of the EM field. Because of this asymmetry, it is possible to construct a genuine matter theory but it is 
difficult to conceive a genuine field theory. 

A most natural derivation of the Casimir force seems to be to calculate the zero-point energy first and then take its 
derivative. The original derivation by Casimir follows this approach. In this type of approach, both field and matter 
contribute to the zero-point energy as in our theory. Nevertheless, one encounters descriptions in the literature to 
the effect that the; Casimir force; is the; hallmark of zero-point energies of the genuine EM field. Perhaps the main 
reason for this (strictly speaking incorrect) recognition lies in the fact that while the system has the region where 
only the EM field (including the vacuum) exists, there is no region where only the matter exists. This also reflects 
the above-mentioned asymmetry between the field and the matter. 
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